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First-principles calculations using density functional theory based on norm-conserving pseudopo- 
tentials have been performed to investigate the Mg adsorption on the Si(001) surface for 1/4, 1/2 
and 1 monolayer coverages. For both 1/4 and 1/2 ML coverages it has been found that the most 
favorable site for the Mg adsorption is the cave site between two dimer rows consistent with the 
recent experiments. For the 1 ML coverage we have found that the most preferable configuration is 
when both Mg atoms on 2x1 reconstruction occupy the two shallow sites. We have found that the 
minimum energy configurations for 1/4 ML coverage is a 2x2 reconstruction while for the 1/2 and 
1 ML coverages they are 2x1. 

PACS numbers: 68.43.Bc, 68.43.Fg 



I. INTRODUCTION 

Adsorption of overlayers of Mg atoms on silicon sur- 
faces has been a subject of growing interest during 
the last decade because of its importance in techno- 
logical applications such as the efficient photocathodes 
and thermionic energy converters. Even though sev- 
eral experimental and theoretical studies have been car- 
ried out for studying the Mg/Si(lll) systems^ only few 
have been reported for the adsorption of Mg on Si(001) 
surfaces The understanding of Mg adsorption 
on Si(001) has become especially an interesting problem 
to investigate because of the possible use of Si substrates 
for the growth of Mg2Si films. More recently, after the re- 
port of superconductivity in MgB 2 , its growth on Si(001) 
surfaces&iii has become a current issue. 

One of the earliest experiments of Mg adsorption on 
Si(001) were carried out by Kawashima et al.^ who have 
obtained LEED and AES results starting from full cover- 
age to lower coverages by allowing for thermal desorption. 
They observed that as the coverage decreased the various 
structural phases appeared in the order lxl, 2x3, 2x2, 
another 2x3, and finally the clean 2x1 which correspond 
to the coverages of 1, 1/3, 1/4, 1/6, and ML, respec- 
tively. They concluded that Mg atoms were adsorbed on 
hollow site (valley-bridge site, in their notation) for each 
coverage mentioned above. 

Hutchison et al., 4 in their STM results for the low 
coverage-Mg/Si(001) case, have reported three types of 
adsorption geometries at room temperature. Type I is 
the most favorable phase and it refers to a single Mg 
atom adsorbed on a cave site. The least frequently ob- 
served type II, being a localized phase, corresponds to 
two Mg atoms adsorbed on cave sites to relieve the stress 
on the chain of type I phases. Type III phase, which 
may occur everywhere but next to a type II phase, con- 
tains multiple Mg adatoms (three at most) adsorbed on 
a combination of adjacent sites. 
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Similar observations also with the use of STM have 
been reported by Kubo et al. 5 for Mg adsorption at room 
temperature. However, they speculated another interpre- 
tation differing from that of Hutchison et al., such that 
according to Kubo et al., type I may correspond to two 
Mg atoms adsorbed on adjacent hollow sites and type 
III may correspond to a single Mg atom adsorbed on a 
pedestal site. In addition, for the annealed case at high 
temperature, they have tentatively suggested two models 
possessing the 2x2 symmetry by incorporating an extra 
Si on the shallow site (see Fig 1) bonding to Mg adatom. 
The first model consists of a single Mg atom adsorbed on 
the neighboring shallow site bonding across the hallow 
site, and second one consists of two Mg atoms adsorbed 
on cave sites on each side of the extra Si atom. 

Kim et al£ have reported in their LEED data that the 
2x2 reconstruction corresponding to 1/2 ML coverage 
occurs at 280° C, whereas for 1/3 ML coverage with 2x3 
symmetry occurs at 390° C showing the dependence of 
surface reconstruction on coverage. Kubo et ali£ did not 
reject the possibility of 1/2 ML as saturation coverage 
for 2x2 surface, either. Cho et ali have performed high 
resolution core-level photoelectron spectroscopy for 2x2 
and 2x3 structures. They did not find any difference 
in the behavior of Mg adsorption on Si(001) surface for 
these two configurations. They have found that the most 
preferable adsorption site for both of these low coverages 
is the cave site between two dimers (bridge site, in their 
notation). 

On the theory side, Khoo and Ong 8 using a semiem- 
pirical self-consistent molecular orbital method have per- 
formed CNDO (complete neglect of differential overlap) 
calculations to investigate the adsorption of Mg atom on 
Si(001) at 1/2 ML coverage. According to their results 
the Mg atom resides on the bridge site above the mid- 
point of Si dimer. 

Wang et al.jii have done first principle calculations 
for the adsorption of another group II element Ba on 
Si(001) at a rather low coverage of 1/16 ML, and they 
have shown that the minimum energy site is the hollow 
site (cave-bridge site, in their notation) where the Ba 
atom is surrounded by four buckled dimers leading to a 
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solitonlike defect in the original c(4 x 2) configuration. 

The aim of this paper is to perform first-principle cal- 
culations for different coverages, i.e., 1/4, 1/2 and 1 ML, 
of Mg on Si(OOl) to get a clear understanding of the ad- 
sorption mechanisms and to investigate the atomic struc- 
ture of the surface covered with magnesium. To the best 
of our knowledge this is the first detailed work to inves- 
tigate theoretically the Mg/Si(001) system for different 
coverages from first principles. 

II. METHOD 

We used pseudopotential method based on density 
functional theory in the local density approximation. 
The self consistent norm conserving pseudopotentials are 
generated by using the Hammann scheme^ which is in- 
cluded in the fhi98PP packaged Plane waves are used 
as a basis set for the electronic wave functions. In order 
to solve the Kohn-Sham equations, conjugate gradients 
minimization method^ is employed as implemented by 
the ABINIT codeM^ The exchange-correlation effects 
are taken into account within the Perdew-Wang schemed 
as parameterized by Ceperly and Alder4£ 

The unit cell included an atomic slab with 8 layers of 
Si plus a vacuum region equal to about 9 A in thick- 
ness. Single-particle wave functions were expanded us- 
ing a plane wave basis up to a kinetic energy cut-off 
equal to 16 Ry. The integration in the Brillioun zone 
was performed using 8 special /c-points sampled using 
Monkhorst-PackiS scheme. Although a couple of cases 
were repeated with 18 special fc-points, no significant im- 
provement has been observed. 

We have used our theoretical equilibrium lattice con- 
stant for the bulk Si (5.405 A) in the surface calculations. 
The bulk modulus for Si is found to be 96 Mb in rather 
good agreement with the experimental result of 97 Mb. 
Our results for Mg bulk lattice parameter is 3.12 A, and 
the c/a ratio is 0.607 that are very close to the experi- 
mental results of 3.20 A and c/a of 0.614, respectively. 
As another check for our pseudopotentials, the lattice 
parameter for Mg2Si was calculated to be 6.28 A which 
is also in close agreement with the experimental value 
of 6.39 AM The energy bands of this compound gives a 
semiconducting gap of E g equal to 0.31 eV. These tests 
suggest that the inclusion of the nonlinear core correc- 
tions to Mg pseudopotential is not needed. 

We have used the 2x2 surface unit cell in our calcu- 
lations to study the adsorption of Mg at a low coverage 
such as 1/4 ML and to include various combinations of 
adsorption sites for the half- and full- monolayer cover- 
ages. The first step was to optimize the clean Si(001)- 
2x2 surface while keeping the two lowest substrate layers 
(out of 8) frozen into their bulk positions and all the re- 
maining substrate atoms were allowed to relax into their 
minimum energy positions. The p(2x2) was found to 
have the lowest total energy with a dimer length of 2.32 
A and a tilt angle of 19° which is in good agreement with 



the experimental dimer-length value of 2.40 ± 0.10 A. 



III. RESULTS AND DISCUSSION 

We have studied the adsorption of Mg atom on the 
Si(001) surface for 1/4, 1/2 and 1 ML starting with the 
reconstructed p(2x2) surface unit cell. We have chosen 
five different sites for adsorption, namely, cave, hallow, 
pedestal, bridge and shallow. The cave site (c) is lo- 
cated above the fourth layer Si, hallow (h) and pedestal 
(p) sites are above the third layer Si, shallow site (s) 
above the second layer Si between two Si dimers, and the 
bridge site (b) is located above the dimer as indicated in 
Fig. 1. The adatom and the upper three monolayers of 
the substrate were then taken to their minimum energy 
configurations by performing structural optimization us- 
ing the Broyden-Fletcher-Goldfarb-Shanno method^ un- 
til the force on each atom reduces to a value less than 25 
meV/A. 
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FIG. 1: Schematics of the adsorption sites for Mg atom on 
Si(001) surface from three different views. The symbols stand 
for: 6=bridge, p=pedestal, Shallow, s=shallow and c=cave. 
(The dimers are shown symmetric here for visual convenience) 

In Table [3 we introduce the structural parameters and 
bond lengths for different adsorption sites as well as the 
adsorption energy E a d (negative of the binding energy of 
the adatom) given by 

Ead = -Bsi(OOl) + "-^Mg - ^Mg/Si(001) (1) 

where -EMg/Si(ooi) is the total energy of the adsorbed sur- 
face, -Esi(ooi) t ne total energy of clean surface, n is the 
number of Mg adatoms in the surface unit cell and Em s 
is the total energy of a single Mg atom with spin polar- 
ization obtained in a separate ab initio calculation using 
the same pseudopotential and the same energy cut off in 
a larger unit cell with a size of about 15 A. 
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A. 1/4 ML Coverage 




(a) (b) 

FIG. 2: Mg overlayer adsorbed on Si(OOl) surface (a) for 1/4 
ML coverage on c site and (b) for 1/2 ML coverage on c-c site 



The adsorption of Mg atom on cave site was found to 
be the most favorable one among the five different ad- 
sorption sites, shown in Fig. 1. The other cases, namely, 
the s site and the b site occupations by an adatom are 
also stable. However, in the case of adsorption of Mg on 
p or h sites, the adatoms reside there temporarily (for 
some optimization steps) before it becomes unstable and 
then they migrate to different sites with lower energies. 
In fact, the adatom on p site migrates to the s site, while 
the adatom on h site migrates to the neighboring c site. 

The adsorption of Mg on c site corresponds to an ad- 
sorption energy of 2.26 eV and s and b sites were found 
to be 0.25 eV and 0.68 eV lower than that of c case, re- 
spectively. A natural pathway for surface diffusion might 
be b — > s — > c. Listed in Table [I] is also the values for 
the metastable case of h site, since it takes longer for the 
h — > c migration, in contrary to the p site with much 
faster p — * s migration since the adatom there is in a 
very unstable environment. 

For all stable cases in 1/4 ML coverage we have found 
that the surface dimers get partially symmetrized. Only 
one dimer ("dimeri" to which Mg adatom bonds) be- 
comes symmetric (with a tilt angle less than ~ 6°), while 
the other one (dimer2 with no charge transfer from Mg 
adatom) is still tilted (e.g., Fig. 2a). The tilting angles 
for dimer 2 in c, b and s cases being 16.0°, 17.1° and 13.0°, 
respectively, are still close to the clean surface value of 
19° while dimeri for all of them are almost flat as seen 
from Table H| 

The height of Mg atom on cave site is found to be 
0.75 A which is in good agreement with the STM re- 
sults of Hutchison et ah 4 Fig. 2a shows for the 1/4 ML 
coverage that the two Si atoms bonded to the same Mg 
adatom are still in an asymmetric environment. The one 
whose neighbors are tilted down is 2.68 A away from Mg 
atom whereas the other has a bondlength of 2.65 A. In 
a different calculation for Mg2Si we have found Mg-Si 
bondlength to be 2.70 A which is close to the back-bond 
lengths of Mg in this environment. 



B. 1/2 ML Coverage 

For the case of 1/2 ML coverage, (i.e., 2 Mg atoms per 
2x2 unit cell) , we have considered the adsorption at com- 
binations of pairs of the five different sites resulting in 15 
different cases among which the most favorable case was 
found to be the c-c sites in our calculations. The struc- 
ture then possesses 2x1 reconstruction with symmetric 
dimers stretched to 2.54 A, shown in Fig. 2b. 

Even though s was stable for the low 1/4 ML coverage, 
we found that s-s is not always stable for 1/2 ML cover- 
age. Out of four s sites in a 2x2 surface unit cell there 
arc three incquivalent s-s geometries as shown in Fig. 3. 
For instance, it is unstable if the other adatom occupies 
the next s site along the [110], i.e., dimer row direction, 
as shown in Fig. 3c. However, if the second adatom were 
added on the other s site along [110], i.e., the direction 
of dimers, as seen in Fig. 3a, this 2x2 configuration will 
be stable and it will have an adsorption energy almost 
equal to that of the c-c case, differing by only 0.01 eV 
per adatom. The last case s-s', with the second s site 
is along the [010] direction as shown in Fig. 3b, is also 
stable with an adsorption energy of 2.14 eV per adatom. 

Although neither of p and h sites are stable at 1/4 



(a) 



(b) 



(c) 




FIG. 3: Adsorption of two Mg atoms on s sites for half cov- 
erage: (a) s-s stable, (b) s-s' stable and (c) s-s" unstable 
configurations 
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TABLE I: The structural parameters (in A) for most of the stable adsorption sites (except for the starred case which is 
metastable) for three different coverages O. The heights of the adatom with respect to the dimers (dx in A) and the adsorption 
energies (E a( j in eV) for these adsorption cases are also presented. The quantities in parenthesis are tilt angles (in degrees) of 
the corresponding dimers. 
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2x2 


2.26 




s 


2.46 (5.5) 2.35 (13.0) 
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2.79 
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2.87 


0.53 


(2)xl 


1.91 
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2.50 


2.65 




2.76 


0.60 1.56 


2x1 


1.74 




c-c-b-b 


2.58 


2.59 




2.69 


0.59 2.36 


2x1 


1.59 




p-p-b-b 


2.56 


2.65 




3.59 


1.31 3.35 


2x1 


1.46 



ML coverage, combination of them at 1/2 ML coverage, 
i.e., the p-h case with one atom occupying a pedestal 
site while the other occupying the neighboring h site was 
found to be stable and it is more favorable than the oc- 
cupation oip-p sites that is also stable. In the case of p-p 
adsorption the adsorbed atoms are four-fold coordinated 
with the neighboring Si dimer atoms with a bond length 
of 2.68 A and with symmetrical dimers stretched to a 
length of 2.52 A. The case of adsorption on b-b was also 
found to be stable with symmetrical dimers but with an 
adsorption energy less than the c-c case by 0.68 eV per 
adatom. 

The other combinations were found to be not stable 
and the adatom migrates to any one of the above stable 
configuration sites. 

In all of the above stable configurations we have found 
that the surface dimers become symmetric (e.g., Fig. 2b) 
in contrast to the 1 /4 ML adsorption case discussed be- 
fore. The dimer lengths for all the stable adsorption sites 
are in the range 2.43-2.55 A except for the case of p-h 



adsorption where the two dimers were stretched to rather 
long values of 2.87 and 2.91 A, respectively. This asym- 
metry makes it slightly a 2x2 rather than a 2x1 like 
most other cases in 1 /2 ML coverage with the exception 
of the s-s case which is naturally a 2x2. As to the Mg- 
Si backbonds, they are in the range of 2.53-2.93 A. In 
Table U the dj^g-Si values listed for the p-h case are for 
the Mg atom on p and h sites bonded to dimer Si atoms, 
respectively. Similarly, the dx values correspond to the 
heights measured from Si dimers of Mg adatoms on p and 
h sites, respectively. 

In addition to the pathway b — > s — > c, we have also 
b p — ► s <-» c pathway in 1/2 ML coverage, where 
the backward migrations in energy take place among sites 
in the trough. Another rule brought by the double occu- 
pancy in this coverage is that the adatoms form chains 
along the dimer rows or along the dimer bonds. The s- 
s'-case looks like an exception to this, however, it is not 
considered to be a chain to first order. Any possible un- 
stable combination a-(3 phase will diffuse into one of the 
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stable phases listed in Table [I] by the adatom on either 
one of a or sites migrating to a new site according to 
the above rules. In the case of b-c phase, since each Si 
atom of all dimers are bonded to two Mg atoms it is not 
stable and consequently, both Mg atoms at b and c sites 
migrate to their corresponding neighbors p and h sites, 
respectively, to reach the stable p-h phase. We should 
emphasize here that this conclusion is not only based on 
the total energy calculations but also on the comparisons 
of the Hellman-Feynman forces and the adatom behavior 
in the other less stable cases. 



For the full coverage adsorption, we have found six 
most probable configurations that are also consistent 
with our findings for the half coverage case. They are 
listed in Table[I]m the order of their adsorption energies. 

The adsorption on s-s-s-s was found to be the most 
favorable one with an adsorption energy of 1.98 eV per 
adatom. In this model we have four Mg atoms occupying 
all four s sites forming Mg lines that run along the dimer 
rows. The surface symmetry of this structure is 2x 1 with 
dimers shortened to 2.38 A as compared to 2.43 A of the 
s-s half coverage case. The Mg-Si backbond of length 
2.73 A compares well with the Mg-Si bond in Mg2Si. 

The next stable case that has less adsorption energy 
(by only 0.05 eV per adatom) is the c-c-s-s in which 
one has two Mg adatoms on s-s along the [110] direc- 
tion, while the other two Mg adatoms occupy the c-c 
sites along the [110] direction as seen in Fig. 4b. In this 
case the reconstructed surface structure possesses a 2 x 2 
symmetry with dimers of length 2.45 A. The Mg-Si back- 
bonds being in the range of 2. 75-2. 79 A also compare well 
with the Mg-Si bonds in silicates. 

The p-p-h-h combination, in which every hallow and 
pedestal sites are occupied, was also found to be stable 
with an adsorption energy differing from the s-s-s-s case 
by 0.07 eV per adatom. This phase is very interesting 
in its symmetry properties. Upon relaxation the dimers 
are removed and it approached to an almost lxl phase 
where p and h become the same. The deviation from 
lxl is very small and is due to a slight shift (0.15 A) of 
only the ex-dimer- member Si atoms in the [110] direction 
opposite to one another, so that the Mg adatom will only 
bond to two Si atoms instead of four. This symmetry 
breaking due to charge transfer from Mg to surface Si 
atoms causes the zigzag bond picture shown in Fig. 4c. 

Next combination is the c-c-p-p case where the ad- 
sorbed atom in p site is four-fold coordinated with the 
neighboring Si dimer atoms with a bond length of 2.65 
A while the other adatom in cave site is two-fold coor- 
dinated with a bond length of 2.76 A giving rise to an 
adsorption energy of 1.74 eV per adatom. This structure 
possesses 2x1 symmetry with a dimer length of 2.50 A. 

The combination c-c-b-b was also found stable with 
adsorption energy equal to 1.59 eV per adatom. In this 




FIG. 4: Three views of the minimum energy structures (after 
relaxation) for (a) Mg/Si(001)-(2x l)-s-s-s-s, (b) Mg/Si(001)- 
(2x2)-c-c-s-s and (c) Mg/Si(001)-(2x l)-p-p-h-h cases for 1 
ML coverage 



C. 1 ML Coverage 
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model we have Mg atoms occupying all the bridge and 
cave sites, forming Mg lines that run orthogonal to the 
dimer rows and resulting in a 2 x 1 symmetry. The dimer 
length was found to be 2.58 A. The difference in heights 
of the Mg adatoms on c and b sites measured from the 
silicon dimers corresponding to 1.77 A undulates the Mg 
chains along the dimer direction. 

The configuration p-p-b-b can be described as double 
half monolayer where two Mg atoms occupy two pedestal 
sites (d^=1.31 A) while the other two Mg adatoms rise 
further up to a height of 3.35 A measured from the center 
of the dimer. This structure possesses 2x1 symmetry 
with a dimer length of 2.56 A and can be described as 
undulated Mg wires located on the Si dimer rows, similar 
but orthogonal to the c-c-b-b case. 

IV. SUMMARY AND CONCLUSION 

We have performed an ab initio total energy calcula- 
tion and geometry optimization for a clean Si(OOl) sur- 
face and that with Mg overlayer of different coverages on 
it. For the 1/4 ML coverage we have found two favorable 
adsorption sites, namely c and s with a relative adsorp- 
tion energy difference of only 0.25 eV per adatom. Our 
results agree well with the room-temperature results of 
Hutchison et al4i and of Cho et ali where they have sug- 
gested that the most favorable site is the cave site. Kubo 
et al£ have suggested for the higher temperature adsorp- 
tion that the s site is favorable, however, they assumed 
an extra Si on the surface bonding to the adatom. Even 
though Kubo et al£ suggested the possibility of single 
adatom occupation on a pedestal site at room temper- 
ature, we did not find any evidence to confirm it. Our 
calculations show that this site is unstable for a single 
adatom adsorption. Occupation of pedestal site has also 
been ruled out by Cho et ali for low coverage. But we 
find that the argument of Hutchison et ali^ that this site 
can be occupied by multiple atoms is in agreement with 
our stable p-p case for 1/2 ML. 

Contrary to the theoretical results^ for Ba adsorption 
on Si(001), we found that the Mg adsorption on h site 
is metastable, i.e., the adatom resides there temporarily, 
altering the dimers into symmetrical ones. However, the 
total energy of this case being higher than the adsorption 
on c site by 0.85 eV per adatom, eventually, it leaves this 
saddle point and migrates to the neighboring cave site. 
We believe the reason for this disagreement with the Ba 
case might be due to the difference in their atomic sizes. 

For the case of half monolayer coverage, our results give 
2x1 reconstruction with the Mg adatoms on c-c sites. 
The high temperature result of Kubo et al.^ for the same 
coverage suggested that it contains two Mg adatoms lo- 
cated nearly on c-c sites but distorted by the existence of 
an extra Si atom which causes the symmetry to be a 2 x 2. 
On the other hand, their 1/2 ML result for low temper- 



ature can be described as two Mg adatoms confined in 
every other cave site and separated towards ft-saddle to 
make it a 2x2. This corresponds to a distorted version 
of our h-h case. Starting from both of these geometries 
we have found them unstable ending up in c-c case with 
a 2 x 1 symmetry. Even though we have started with a 
2x2 unit cell for each case in our calculations we have 
ended up with 2x1 reconstruction for the 1/2 ML cov- 
erage except for the s-s cases, which are naturally 2x2, 
and the p-h case which is slightly a 2x2 distorted from 
2x1. This slight symmetry breaking shows up itself in 1 
ML coverage case of p-p-h-h as well. 

The bridge site b-b was found to be stable in accor- 
dance with the theoretical result of Khoo and OngS for 
this coverage, even though it is the least favorable case. 

For the case of full monolayer coverage, we have found 
that the adsorption model in which all the Mg atoms oc- 
cupy all the s sites with 2x1 construction is more favor- 
able than 2x2 reconstruction in which Mg atoms occupy 
two shallow and two cave sites having an adsorption en- 
ergy difference of only 0.05 eV per adatom. The next 
favorable case having 0.07 eV per adatom smaller ad- 
sorption energy than the most preferable s-s-s-s case is 
the p-p-h-h configuration with a 2x1 symmetry, but just 
slightly away from lxl geometry as seen in Fig. 4c. This 
may suggest that more than 1 ML coverage might lead 
to a 1 x 1 phase consistent with the LEED results of Kim 
et al£ who reported that the lxl phase corresponds to 
2 ML coverage. 

The model of adsorption that can be concluded from 
our theoretical results is that for the low coverage, Mg 
atoms will be adsorbed mostly on c sites and less likely on 
s sites, which continues until the coverage reaches to 1/4 
ML. As the coverage exceeds 1/4 ML, the Mg atoms are 
adsorbed on the c-c sites and on the s-s sites along [110] 
with almost equal probabilities until the surface becomes 
saturated at 1/2 ML, where it is going to have domains 
of c-c and s-s phases with 2x1 and 2x2 reconstructions, 
respectively. Further increase in coverage will result into 
s-s-s-s, c-c-s-s and even p-p-h-h domains with 2x1, 2x2 
and almost lxl reconstructions, respectively. 

Using ab initio total energy calculations we have per- 
formed detailed investigation of the atomic geometry of 
Mg adsorbed Si(001) surface for 1/4, 1/2 and full mono- 
layer coverages. We have investigated the adsorption 
sites, the energetics and the reconstruction of Si(001) sur- 
face after Mg adsorption. We have found that our results 
agree well with the existing experimental observations. 
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